505 


ACES  JOURNAL,  Vol.  29,  No.  7,  JULY  2014 


Electromagnetic  Scattering  by  Multiple  Cavities  Embedded  in  the 

Infinite  2D  Ground  Plane 


Peijun  Li  ^  and  Aihua  W.  Wood  ^ 

^  Department  of  Mathematies 
Purdue  University,  West  Lafayette,  fN  47907,  USA 
lipeijun@math.purdue .  edu 

^Department  of  Mathematies  and  Statisties 
Air  Foree  Institute  of  Teehnology,  WPAFB,  Dayton,  OH  45433,  USA 
aihua.wood@afit.edu 


Abstract  —This  paper  is  eoneerned  with  the 
mathematieal  analysis  and  numerieal  eomputation 
of  the  eleetromagnetie  wave  seattering  by  multiple 
open  eavities,  whieh  are  embedded  in  an  infinite 
two-dimensional  ground  plane.  By  introdueing  a 
new  transparent  boundary  eondition  on  the  eavity 
apertures,  the  seattering  problem  is  redueed  to  a 
boundary  value  problem  on  the  two-dimensional 
Helmholtz  equation  imposed  in  the  separated 
interior  domains  of  the  eavities.  The  existenee  and 
uniqueness  of  the  weak  solution  for  the  model 
problem  is  studied  by  using  a  variational  approaeh. 
A  bloek  Gauss-Seidel  iterative  method  is 
introdueed  to  solve  the  eoupled  system.  Numerieal 
examples  are  presented  to  show  the  effieieney  and 
aeeuraey  of  the  proposed  method. 

Index  Terms  -  Eleetromagnetie  eavity,  finite 
element  method,  Helmholtz  equation,  variational 
formulation. 

I.  INTRODUCTION 

A  eavity  is  referred  to  as  a  loeal  perturbation 
of  the  infinite  ground  plane.  Given  the  eavity 
strueture  and  an  ineident  wave,  the  seattering 
problem  is  to  prediet  the  eleetromagnetie  field 
seattered  by  the  eavity.  It  has  been  extensively 
examined  by  researehers  for  the  time-harmonie 
analysis  of  eavity-baeked  apertures  with 
penetrable  material  filling  the  eavity  interior  [14- 
16,  18,  28].  Mathematieal  analysis  of  the  problem 
ineluding  overfilled  eavities,  where  the  aperture  is 
not  planar  and  may  protrude  the  ground  plane,  ean 


be  found  in  [1-4,  17,  19-24,  27].  A  lot  of  work  has 
been  devoted  to  solve  the  problem  by  various 
numerieal  methods  ineluding  finite  element,  finite 
differenee,  boundary  element,  and  hybrid  methods 
[5,  7,  8,  11,  12,  25,  26,  29,  30].  All  the  model 
problems  have  been  foeused  on  a  single  eavity, 
whieh  may  limit  the  applieation  of  the  problem  in 
industry  and  military.  This  paper  aims  to  extend 
the  single  eavity  model  to  a  more  general  multiple 
eavity  model,  and  analyze  and  develop  numerieal 
methods  for  the  assoeiated  boundary  value 
problem. 

In  this  paper  we  foeus  on  the  Transverse 
Magnetie  polarization  (TM),  where  the  modeling 
equation  is  the  two-dimensional  Helmholtz 
equation.  Based  on  Fourier  transform,  a  nonloeal 
transparent  eondition  is  introdueed  on  the  aperture, 
whieh  eonneets  the  eleetrie  field  in  eaeh  individual 
eavity.  By  using  the  boundary  eondition,  we 
reduee  the  multiple  eavity  problem  into  a 
boundary  value  problem  imposed  in  the  interiors 
of  the  eavities.  The  existenee  and  uniqueness  of 
the  weak  solution  for  the  model  problem  is  studied 
by  using  a  variational  approaeh.  A  bloek  Gauss- 
Seidel  iterative  method  is  introdueed  to  solve  the 
eoupled  system,  where  only  a  single  eavity 
problem  needs  to  be  solved  at  eaeh  iteration.  Thus, 
it  is  applieable  of  any  effieient  single  eavity  solver 
to  the  multiple  eavity  problem.  Numerieal 
examples  are  presented  to  show  the  effieieney  and 
aeeuraey  of  the  proposed  method.  We  refer  to  [9, 
10,  13]  for  numerieal  methods  to  solve  a  related 
multiple  obstaele  seattering  problem. 
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The  paper  is  outlined  as  follows.  In  Seetion  2, 
a  mathematieal  model  for  the  single  eavity 
problem  is  introdueed,  the  variational  formulation 
is  presented,  and  the  uniqueness  and  existenee  of 
the  solution  are  examined.  Seetion  3  is  devoted  to 
the  study  of  the  solution  for  the  multiple  eavity 
problem.  The  major  new  ingredient  is  the 
introduetion  of  a  transparent  boundary  eondition. 
Seetion  4  addresses  the  numerieal  implementation 
and  examples  are  shown  to  illustrate  the  method. 
The  paper  is  eoneluded  with  some  general  remarks 
and  direetions  for  future  researeh  in  Seetion  5. 

II.  SINGLE  CAVITY  SCATTERING 

In  this  seetion,  we  study  a  single  eavity 
problem,  whieh  is  intended  to  serve  as  a  basis  for 
the  multiple  eavity  problem. 

A.  Model  problem 

We  foeus  on  a  two-dimensional  geometry.  The 
medium  is  assumed  to  be  non-magnetie  and  has  a 
eonstant  magnetie  permeability;  i.e.,  ju  =  where 

ju^  is  the  magnetie  permeability  of  vaeuum.  The 

medium  is  eharaeterized  by  the  dieleetrie 
permittivity  s . 

As  shown  in  Fig.  1,  an  open  eavity  Q 
enelosed  by  the  aperture  r  and  the  wall  S  ,  is 
plaeed  on  a  perfeetly  eondueting  ground  plane  F"^ . 
Above  the  flat  surfaee  {y  =  0}  =  F  u  F%  the  medium 
is  assumed  to  be  homogeneous  with  a  positive 
dieleetrie  permittivity  .  The  medium  inside  the 

eavity  is  inhomogeneous  with  a  variable  dieleetrie 
permittivity  £{x,y) .  Assume  further  that: 

s  e  Re  ^  >  0,  Im^  >  0 . 

For  the  TM  polarization,  the  magnetie  field  is 
transverse  to  the  invariant  direetion.  The  time- 
harmonie  Maxwell  equations  ean  be  redueed  to  the 
two-dimensional  Helmholtz  equation: 

At/ +  =  0  in  Q  u  .  (1) 

The  total  field  satisfies  the  boundary  eondition: 

u  =  0onrKjS,  (2) 

where  is  the  wavenumber  and  co  is  the 

angular  frequeney. 

Let  an  ineoming  plane  wave  u' 
be  ineident  on  the  eavity  from  above,  where  0  is 
the  ineident  angle  with  respeet  to  the  positive 
y  axis,  and  wavenumber  of  the 

free  spaee. 


Denote  the  referenee  field  t/"^^  as  the  solution 
of  the  homogeneous  Helmholtz  equation  in  the 
upper  half  spaee: 

=  0  in  Rl,  (3) 

together  with  boundary  eondition: 

t/^^^  =  0on  r  uF.  (4) 

It  ean  be  shown  from  (3)  and  (4)  that  the  referenee 
field  eonsists  of  the  ineident  field  and  the  refleeted 
field: 


where  ^ 

The  total  field  is  eomposed  of  the  referenee 
field  and  the  seattered  field: 

u  =  +  u\ 

It  ean  be  verified  from  (1)  and  (3)  that  the 
seattered  field  satisfies: 

=0m  Rl.  (5) 

In  addition,  the  seattered  field  is  required  to  satisfy 
the  radiation  eondition: 


p^oc  Qp 

To  deseribe  the  boundary  value  problem,  we 
need  to  introduee  some  funetional  spaees.  For 
//GZ^(r'^uF),  whieh  is  identified  with  l}{R),  we 
denote  by  u  the  Fourier  transform  of  u  defined  as: 


u{^)  =  j  u{x)e'^^dx. 

R 

Using  Fourier  modes,  the  norm  on  the  spaee 
l}{R)  ean  be  eharaeterized  by: 


1  1 


^\uf  dx 

2 

\\ufd^  ' 

_R 

_R 

Denote  the  Sobolev  spaee: 

=  Z"  (Q),  for  1 5  |<  1}  , 

and  the  traee  funetional  spaee: 

H^(R)  =  {«  e  e{R)  + 

R 


whose  norm  is  defined  by: 


\u\\ 


'H\R) 


\{\+er\ufd^ 

R 


By  taking  the  Fourier  transform  of  (5)  with 
respeet  tox ,  we  obtain: 

Since  the  solution  of  (7)  satisfies  the  radiation 
condition  (6),  we  deduce  that; 


=  (8) 
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where 


m)  =  \ 


{Ki-eyiov\^^K„ 

i{e-K:lYior\^\>K,. 


Taking  the  inverse  Fourier  transform  of  (8),  we 
find  that: 


R 

Taking  the  normal  derivative  on  V  uT,  whieh 
is  the  partial  derivative  with  respeet  to  y  ,  and 
evaluating  at  3;  =  o  yield: 

dX(x,y)  1^.0=  I  ij3i^)u(4,0)e-‘^d4.  (9) 

R 


For  given  u  on  r'^uF  ,  define  the  boundary 
operator  T : 

Tu  =  \  iP{^)u{^,0)e-‘^d^,  (  i  9) 

R 

whieh  leads  to  the  transparent  boundary  eondition 
for  the  seattered  field  on  F^  u  F : 

dy{u-u''^)  =  T(u-u^^^. 

Equivalently  we  have  a  transparent  boundary 
eondition  for  the  total  field: 

=  7h  +  gOn  F""  uF ,  (11) 

where 

g  =  =  -liK^ 

It  ean  be  shown  that  the  boundary  operator  is 

1  3 

eontinuous  from  H^{K)Xo  H~^{R).  Furthermore,  it 
has  the  following  properties  whieh  ensures  the 
uniqueness  of  the  solution  of  the  single  eavity 
problem. 


Lemma  1.  Let  u  e  H^{R) .  It  holds  that  Rq(Tu,u)  <  0 
and  lm(Tu,u)>0  •  Furthermore,  if  u  is  an 


analytical  function  with  respect  to^,  Re(rw,w)  =  0 
or  \m{Tu,u)  =  0  implies  u  =  0. 


To  derive  a  transparent  boundary  eondition  for 
the  total  field  on  the  aperture  F,  we  need  to  make 
the  zero  extension  as  follows:  for  any  given  u  on 
F,  define 


u(x)  = 


[  W,X  G  F, 
|o,xgF^ 


The  zero  extension  is  eonsistent  with  the  problem 
sinee  the  ground  plane  is  a  perfeetly  eleetrieal 
eonduetor.  Based  on  the  extension  and  the 
transparent  boundary  eondition  (11),  we  have  the 


transparent  boundary  eondition  for  the  total  field 
on  the  aperture: 

dyU  =  Tu+gOnr.  (12) 

1"  F  F" 


S 

Fig.  1.  The  problem  geometry  of  a  single  eavity. 

B.  Well-posedness 

Define  a  traee  funetional  spaee: 

H^(T)  =  {u:u  eH^(R)}, 

1 

whose  norm  is  defined  as  the  h^(R)  norm  for  its 
extension;  i.e.. 

Define  a  dual  paring: 

(^,v)r  =\uv. 
r 

This  dual  paring  for  u  and  v  is  equivalent  to  the 
sealar  produet  in  F{R)  for  their  extensions;  i.e., 

_1  1 

Denote  by  //  2(r)the  dual  space  of  H^{T)',  i.e., 

_1  1 

H  =  .  The  norm  on  this  spaee  is 

eharaeterized  by: 

II II  . 

IHL-3,r,=  sup 

Introduee  a  spaee: 

//^(Q)  =  {u^H\a):u  =  0onS}, 
whieh  is  a  Hilbert  spaee  with  the  usual  norm. 

Multiplying  a  test  funetion  v  on  both  sides  of 
(1)  and  using  the  boundary  eonditions  (2)  and  (12), 
we  may  deduee  a  variational  problem:  find  u  sueh 
that 

a{u,v)  =  (^g,v)y  forall  vg//^(Q),  (13) 

where  the  sesquilinear  form  is: 

a{u,  v)  =  I ( Vt/  •  Vv  -  K^uv)  -  (Tu,  v)^ .  ( ^ 

Q 


Theorem  1.  The  variational  problem  (13)  has  a 
unique  weak  solution  in  ond  the  solution 
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satisfies  the  estimate: 

where  C  is  a  positive  constant. 

Proof:  Decompose  the  sesquilinear  form  (14) 
into  a  =  a, where 

a^(u,v)  =  j  Vi/  •  Vv  -  {Tu,v)  ’ 

and 

a^iuw)  -^K^uv. 

Q 

We  conclude  that  from  Lemma  1  and  Poincare 
inequality  that  a^is  coercive  from: 

Rea^(w,w)  =  j|  Vw  p  -Re(7w,v)  >j|  Vw  | 

Q  Q 

foraUMe//‘(Q). 

Next  we  prove  the  compactness  of <22.  Define 
an  operator  k  :  Z"(Q)  ^  7f^(Q)  by: 

{Ku,  v)  =  {u,  v)  for  aU  v  g 
which  explicitly  gives  that  for  all  v  g  hI{Q)  , 

J  VKu  ■  Vv  -  {TKm,  V )  =  J  k\v. 

n  Q 

Using  the  coercivity  of  and  the  Lax-Milgram 
lemma,  it  follows  that: 

NL.(n)^^IHL^(n)- 

Thus,  K  is  bounded  from  P{Q)  to  h\Q)  and 
H\Q)  is  compactly  imbedded  into  Z^(Q).  Hence, 
K :  Z^(Q)  ^  Z^(Q)  is  a  compact  operator. 

Define  a  function  wgL^(Q)  by  requiring 
wsHliO)  and  satisfying: 

a^(w,v)  =  for  aU  v  g  7V^(Q). 

It  follows  from  the  Lax-Milgram  lemma  again  that: 

Using  the  operator  K  ,  we  can  see  that  the 
variational  problem  (13)  is  equivalent  to  find 
ueP{D.)  such  that: 

(I-K)u  =  w.  (17) 

It  follows  from  the  uniqueness  result  and  the 
Fredholm  alternative  that  the  operator  I-K  has  a 
bounded  inverse.  We  then  have  the  estimate: 

Combining  (15)-(17),  we  deduce  that: 


ll%-(n)^INL.(n)+IM 


<  C  b  o  +  \\w\\  , 


which  completes  the  proof 

III.  MULTIPLE  CAVITY  SCATTERING 

As  shown  in  Fig.  2,  we  consider  a  situation  of 
n  cavities,  where  the  multiple  open  cavities 
Qp...,Q^  enclosed  by  the  apertures  L,...,r^  and  the 
walls  are  placed  on  .  Above  the  flat 

surface  {j;  =  0}  =  U u---ur^ uFS  the  medium  is 

assumed  to  be  homogeneous  with  a  positive 
dielectric  permittivity  .  The  medium  inside  the 

cavity  is  inhomogeneous  with  a  variable 
dielectric  permittivity  sfx.y)  ,  which  satisfies 
Re^.  >0,  lm^.^0  for  J  = 

We  consider  the  two-dimensional  Helmholtz 
equation  for  the  total  field: 

Au-\-K:f  =  0,mQ.^  u---uQ^  (1^) 

together  with  the  boundary  condition: 

u  =  0,onS,u--uS^ur.  (19) 

Let  the  plane  wave  u'  be  incident  on  the  cavities 
from  above.  The  total  field  u  is  consisted  of  the 
incident  field  u' ,  the  reflected  field  i/" ,  and  the 
scattered  field  u'  ,  where  the  scattered  field  is 
required  to  satisfy  the  radiation  condition  (6). 

To  reduce  the  problem  into  the  bounded 
domains  =  1,...,^,  we  need  to  derive  a 

transparent  boundary  condition  on  .  Rewrite 
(18)-(19)  into  n  single  cavity  scattering  problem: 

i^Uj+K]u.=0m€ij,  ^20) 

Uj  =  0  on  Sj, 

where  .  If  ^  is  the  solution  of  (18)-(19) 

and  Uj  is  the  solution  of  (20),  respectively,  then  we 
have  for  j  = 

For  u.{xfi)^  define  its  zero  extension: 

fw  (x,0)  forx  G  F., 


w  .(x,0) : 


Oforx  G  R  \  F.. 


For  the  total  field  u ,  define  its  extension: 

m  K(-r,0)forxGF 
u{xSj)  =  < 

\  Oforx  gF^ 

It  follows  from  the  definition  of  the  extensions  that 
we  have: 

n 

u  =  'Yfij  on U  u  •  •  •  u F^  u F"". 

7=1 

Repeating  the  same  steps  as  those  for  the  single 
cavity  problem,  we  have  the  following  transparent 
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boundary  condition  for  the  extended  field: 

a/  =  7w+gonriU---ur„ur,  (21) 
which  gives  the  transparent  boundary  for  Uj  • 

dUj=  Tu.  +  V  +  g  on r  . 

y  j  j  I  ^  j 

i=U^j 

As  we  can  see  from  (22),  the  boundary  condition 
for  UjJ  =  \,...,n  is  coupled  with  each  other,  which 

is  the  major  difference  between  the  single  cavity 
problem  and  the  multiple  cavity  problem. 

Next  we  present  a  variational  formulation  for 
the  multiple  cavity  problem.  Denote 
Q  =  Qj  u  •  •  •  u  ,  r  =  u  •  •  •  u  and  S  = 

Define  a  trace  functional  space: 

Its  norm  is  characterized  by: 

HI  i  =ZhlH  • 

y=i  7?2(r.) 

Denote  // 2(r)  =  // ^^rjx-.-x 

1 

the  dual  space  of  h^{T).  The  norm  on  the  space  is 
characterized  by: 

n  2 

wr  i  =zhii .  • 

Introduce  the  space: 

//^(Q)=/4^(Q0x...x/4(qj, 
which  is  a  Hilbert  space  with  norm  characterized 
by: 

2  "  " 

HH(o)  =  ZK-|| 

7=1  H\n^) 

Similarly,  we  may  obtain  the  variational 
formulation  for  the  multiple  cavity  problem:  find 
u  G  H\{a)  with  u  .=u\^  such  that 


:H  2(r^)  ,  which  is 


a{u,v)  =  foraU  v  g  (^3) 

7=1 

where  the  sesquilinear  form  is: 

a{u,v)  =  ^\(yUj-VVj  - 

7=1  Q^.  7=1  «=1 

We  have  the  following  well-posedness  result. 
The  proof  is  similar  in  nature  as  that  of  the  single 
cavity  problem  and  is  omitted  here  for  brevity. 

Theorem  2.  The  variational  problem  (23)  has  a 
unique  weak  solution  in  and  the  solution  satisfies 
the  estimate: 


ll^ll/7*(Q)  “  2(r)  ’ 

where  C  is  a  positive  constant. 

pc 


pc 


5i 


5. 


Fig.  2.  The  problem  geometry  for  multiple  cavities. 

IV.  NUMERICAL  EXPERIMENTS 

In  this  section,  we  discuss  the  computational 
aspects  and  present  some  examples  for  the 
multiple  cavity  problem. 

A.  Finite  element  formulation 

Let  Mj  be  a  regular  conforming  triangulation 

of  Q  and  V.  (Q  )  be  the  conforming  linear 
finite  element  space  over  Mj  •  Denote 
U  =  Fi  X  •  •  •  X  .  The  finite  element  approximation  to 
the  multiple  cavity  problem  is  to  find  with 
f  G  V.  such  that 


a{u^  ’  ^^  )  “  L  )r  ^ 

7=1  ^ 

where  the  sesquilinear  form 

a(«\v‘)=x|(v«*.vv;-^>;v;)-xx(7«;,v,^). 

7=1  Q^.  7=1  i=l 

For  any  i<  j  <n,  we  denote  by  p.  the  set  of 
vertices  of  Mj,  which  are  not  on  the  cavity  wall  Sj ? 
and  let  (p.{r)^Vj  be  the  nodal  basis  function 
belonging  to  vertex  r^P.  •  Using  the  basis 
functions,  the  solution  of  (24)  is  represented  as: 

The  discrete  problem  (24)  is  equivalent  to  the 
following  system  of  algebraic  equations: 

(25) 

where 


A  = 


AU=G, 

4-^,1 

~  ^2,1 

•  -4, 

-42 

> 

..  1 

KJ 

•  -4.2 

-4,„ 

-4,„  •• 

•  4-4. 
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Hence,  each  Uj  is  an  unknown  vector  whose 

entries  are  Uj{r)  =  u^.{r)  for  all  r^Pj  ,  Aj  is  the 

stiffness  matrix  for  the  discrete  problem  and  its 
entries  are  defined  by: 

Aj  {r,r')=  ^[V(pj{r)-V  (pj  (r' )  -  K^cpj  {r)(p.  (r' )], 


for  all  r,r'(^Pj  ■  The  entries  of  5  .  are  defined  by: 

B.Xr,r')  =  iT^J^)  for  all  r,r’GP.  nH, 
and  the  entries  of  each  vector  g  are  given  by: 

gjiY  =  {g.cPjir))^  for  aU  r  e  Pj  n  Tj. 


A  block  Gauss-Seidel  method  is  adopted  to 
solve  (25).  Given  ,  define  U^’'\k>l  by  the 
solution  of  the  following  system  of  equations: 

7-1 


(k) 


(26) 


i^j+l 

The  block  Gauss-Seidel  iteration  (26)  is  equivalent 
to  apply  the  finite  element  method  to  solve  the 
following  problem:  let  uf^  =  0  ?  define  f^^ 


A:  >  1  by  the  solutions  of  the  decoupled  equations 

/U  . . . 

(27) 


Au^  +  ^  llki  Q^., 

uY  =  0  on  Sj , 


d^uf^  =  +  g  onT., 

/=!  7+1 

for  y  =  1,...,^.  Therefore,  we  only  need  to  solve  a 
single  cavity  problem  (27)  at  each  iteration. 


B.  Transparent  boundary  condition 

The  transparent  boundary  conditions  (11)  and 
(22)  are  not  convenient  to  be  implemented 
numerically.  We  take  an  alternative  and  equivalent 
transparent  boundary  condition  [7]. 

Let 

G(r,  r' )  =  T  [//«  - //« 

be  the  Green  function  of  the  two-dimensional 
Helmholtz  equation  in  the  upper  half  space,  where 
is  the  Hankel  function  of  the  first  kind  with 

order  zero;  r  =  (x,y),A=(x\y),p=\r-A\,p=\r-F\, 
and  p=(x\-y)  is  the  image  of  r'  with  respect  to 


the  real  axis.  By  the  Green’s  theorem  and  the 
radiation  condition,  we  obtain: 

dy (x,0)  =  ^  f  — ^ 77® {K,\r-r'\y (r' fi)dr' , 

2  ^  I  r  -  r'  I 

where  is  the  Hankel  function  of  the  first  kind 
with  order  one.  Hence,  the  alternative  boundary 
condition  is: 

dyU  =  Tu  +  gonT,  (28) 

where  the  boundary  operator  T  is  defined  as: 

Tu=^-^\ — i — HY{Kyr-F\)u{r\0)dr\  (^9) 

Here  the  integral  is  understood  in  the  sense  of 
Hadamard  finite-part.  For  multiple  cavities  with 
apertures  u  •  •  •  u  ,  the  boundary  operator  is 

defined  as: 


7^^  =  ^  Z  j  r^vT  (^0  k  -  l)^(^'  .9)^^’  • 

z  j=i  r  J  ^  ^  I 

The  boundary  operator  (29)  can  be  approximated 
by: 

Tu{Xi) »  Y^gMXk,^\ 

k=\ 

where 


and 


Reg*  =-/,t 

Img,,  = 


^0 1  -  Xj  I 


YMo\Xi-x,  I), 


JMAXj  -Xj  I) 


X,  -  X, 


T(i_]n2)for|/-/t|=  1, 

h, 

2 

=i  - for  I  /  -  ^  1=  0, 

K 

T  In  — — —  for  I  /  -  7  |>  2, 

K  |l-^l  -1 
where  is  the  step  size  of  the  partition  for  the 


cavity  aperture  r,  L  7^  are  Bessel  functions 
of  the  second  and  first  kind  with  order  one, 
respectively.  Therefore,  the  boundary  integral 
{Tu,v)  in  the  weak  formulation  for  the  cavity 

problem  can  be  approximated  by  any  numerical 
quadratures. 


C.  Numerical  examples 

The  physical  parameter  of  interest  is  the  Radar 
Cross  Section  (RCS),  which  is  defined  by: 

C7  =  —  \P{,l>)f- 

«'o 
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Here  ^  is  the  observation  angle  and  P  is  the  far- 
field  eoeffieient  given  by: 

When  the  ineident  and  observation  direetions  are 
the  same,  cj  is  ealled  the  baekseatter  RCS,  whieh 
is  defined  by: 

Baekseatter  RCS  (^)  =  10  log^o  cr(^)dB. 

Example  1.  Consider  a  plane  wave  seattering  from 
a  reetangular  eavity  with  1  meter  wide  and  0.25 
meters  deep  at  normal  ineidenee;  i.e.,  6^  =  0.  Two 
different  eases  are  eonsidered:  an  empty  eavity 
with  k  =  Kq  and  a  eavity  filled  with  a 

homogeneous  medium  with  =  ^*^(4  +  /) .  These 

two  eases  have  been  eonsidered  as  standard  test 
problems  in  [14].  The  Reetangular  domain 
[-0.5, 0.5]  X  [-0.25, 0.0]  is  first  divided  into  160x40 
small  equal  reetangles  and  then  eaeh  small 
reetangle  is  subdivided  into  two  equal  triangles. 
Numerieal  results  are  obtained  by  using  a  linear 
finite  element  over  triangles  at  the  wavenumber 
Kq  =  Ik  •  Figures  3  and  4  show  the  magnitude  and 

the  phase  of  the  total  field  on  the  aperture  at  the 
normal  ineidenee,  the  baekseatter  RCS  for  the 
empty  eavity  and  the  filled  eavity,  respeetively. 
We  observe  the  eoineidenee  of  the  numerieal 
results  obtained  in  [19]  (eireled)  and  our  numerieal 
method  (solid  line). 


Fig.  3.  The  magnitude,  phase,  and  baekseatter 
RCS  of  the  total  field  for  Example  1  of  the  empty 
eavity. 


Fig.  4.  The  magnitude,  phase,  and  baekseatter 
RCS  of  the  total  field  for  Example  1  of  the  filled 
eavity. 

Example  2.  Consider  the  normal  ineidenee  of  a 
plane  wave  onto  two  identieal  reetangular  eavities. 
Eaeh  eavity  is  1  meter  wide  and  0.25  meters  deep; 
they  are  1  meter  distanee  away  from  eaeh  other. 
The  two  reetangular  domains  are  given  as  follows: 
eavity  one:  [-i .5,  -  0.5]  x  [-0.25, 0.0], 
eavity  two:  [0.5, 1.5]  x  [-0.25, 0.0]. 

Eaeh  reetangular  domain  is  divided  into  160x40 
small  equal  reetangles  and  then  eaeh  small 
reetangle  is  subdivided  into  two  equal  triangles. 
Three  types  of  eavities  are  eonsidered:  (type  one) 
two  empty  eavities  with  \  (type  two) 

two  filled  eavities  with  ^^^2^4  +  /) ;  (type 

three)  one  empty  eavity  with  and  one  filled 

eavity  with  =  a*q(4  +  /)-  Figures  5,  6  and  7  show 

the  magnitude  and  the  phase  of  the  total  field  on 
the  apertures  at  the  normal  ineidenee  and  the 
baekseatter  RCS  for  the  type  one,  type  two  and 
type  three  eavities,  respeetively.  These  numerieal 
results  are  obtained  by  the  bloek  Gauss-Seidel 
iterative  method.  To  show  the  eonvergenee  of  the 
iterative  method,  we  define  the  error  between  two 
eonseeutive  approximations: 

^  i<j<k\\  2  J  llL(r.)’ 

where  k  is  the  number  of  iteration.  Figure  8  shows 
the  error  of  two  eonseeutive  approximations 

against  the  number  of  iterations  for  all  three  types 
of  eavities.  It  ean  be  seen  from  Fig.  8,  that  more 


LI,  WOOD:  ELECTROMAGNETIC  SCATTERING  BY  MULTIPLE  CAVITIES 


512 


number  of  iterations  are  needed  for  the  type  one 
eavities  to  reaeh  the  same  level  aeeuraey  as  the 
other  two  types  of  eavities.  The  reason  is  that  the 
eavity  for  either  type  two  or  type  three  is  filled 
with  eomplex  medium,  whieh  aeeounts  for  the 
absorption  of  the  energy,  and  thus,  the  damping  of 
the  amplitude  of  the  field. 


Fig.  5.  The  magnitude,  phase,  and  baekseatterer 
RCS  for  Example  2  of  the  type  one  eavity. 


Fig.  6.  The  magnitude,  phase,  and  baekseatterer 
RCS  for  Example  2  of  the  type  two  eavity. 


Fig.  7.  The  magnitude,  phase,  and  baekseatterer 
RCS  for  Example  2  of  the  type  three  eavity. 


Fig  8.  Convergenee  of  the  Gauss-Seidel  iteration 
for  Example  2. 

Example  3.  Consider  the  seattering  of  a  triple 
eavity  model.  Let  a  plane  wave  be  ineident  onto 
three  identieal  reetangular  eavities  at  the  normal 
direetion.  Eaeh  eavity  is  1  meter  wide  and  0.25 
meters  deep;  there  are  1  meter  distanee  away  from 
eaeh  other.  The  three  reetangular  domains  are 
given  as  follows: 

eavity  one:  [-2.5,  - 1.5]  x  [-0.25, 0.0], 
eavity  two:  [-0.5, 0.5]  x  [-0.25, 0.0], 
eavity  three:  [I.5,  2.5] x [-0.25,  0.0] . 

Again,  eaeh  reetangular  domain  is  divided  into 
160x40  small  equal  reetangles  and  then  eaeh  small 
reetangle  is  subdivided  into  two  equal  triangles. 
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Cavities  one  and  three  are  filled  with  the  same 
homogeneous  medium  with  =Kl(4  +  i)  and 

eavity  two  is  an  empty  eavity  with  v  =  .  Figure  9 

shows  the  magnitude  and  the  phase  of  the  total 
field  on  the  apertures  at  the  normal  ineidenee  and 
the  baekseatter  RCS. 


Fig  9.  The  magnitude,  phase,  and  baekseatter  RCS 
of  the  total  field  for  Example  3. 

V.  CONCLUSION 

We  studied  the  problem  of  eleetromagnetie 
seattering  by  multiple  eavities  embedded  in  the 
infinite  two-dimensional  ground  plane.  The 
seattering  problem  was  redueed  into  a  boundary 
value  problem  by  introdueing  a  transparent 
boundary  eondition.  Based  on  the  variational 
formulation,  we  proved  the  uniqueness  and 
existenee  of  the  weak  solution  for  the  model 
problem.  We  employed  a  bloek  Gauss-Seidel 
iterative  method  to  deeouple  the  eoupled  system 
arising  from  the  multiple  interaetion  among 
eavities.  At  eaeh  step  of  iteration,  it  required  to 
solve  only  a  single  eavity  problem.  Three 
numerieal  examples  were  eonsidered,  a  single 
eavity,  two  eavities  and  three  eavities,  with  and/or 
without  filling.  The  results  show  the  eonvergenee 
of  the  bloek  Gauss-Seidel  iterative  method  for  the 
examples.  We  point  out  some  future  direetions 
along  the  line  of  our  present  work.  The  first  is  to 
analyze  the  eonvergenee  of  the  Gauss-Seidel 
iterative  method  and  investigate  the  parameters, 
sueh  as  separation  distanee  among  eavities. 


wavenumber  and  eavity  size,  whieh  requires 
further  mathematieal  analysis  of  the  stability  of  the 
eavity  seattering  problem  [6].  Another  projeet  is  to 
study  the  multiple  overfilled  eavity  problem  and 
the  model  problem  of  three-dimensional  Maxwell 
equations. 
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